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Abstract 



A unified presentation of the perturbation and variational methods for the 

generalized statistical mechanics based on Tsallis entropy is given here. In 

the case of the variational method, the Bogoliubov inequality is generalized 

in a very natural way following the Feynman proof for the usual statistical 

mechanics. The inequality turns out to be form-invariant with respect to the 

entropic index q. The method is illustrated with a simple example in classical 

mechanics. The formalisms developed here are expected to be useful in the 

discussion of nonextensive systems. 
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Nonextensive effects are common in many branches of the physics, for instance, anoma- 
lous diffusion |[]-§J, astrophysics with long-range (gravitational) interactions [0-fJ], some 
magnetic systems [flO] p~^| , some surface tension questions |T3yi4f . These examples indi- 
cate that the standard statistical mechanics and thermodynamics need some extensions. 
In this direction, a theoretical tool based in a nonextensive entropy (Tsallis entropy) |T5j 
has successfully been applied in several situations, for example, Levy-type anomalous su- 



perdiffusion |16| , Euler turbulence fll7fl , self-gravitating systems [p^7H21| , cosmic background 



radiation f2"2fl ) peculiar velocities in galaxies 23|1 , linear response theory |24] and eletron 



phonon interaction ||25|| , and ferrofluid-like systems ||26|| . Lavenda and coworkers |27j stress 



that any newly proposed entropy [28] must have "concavity" property for it to be correct. 



Tsallis [TjJ and Mendes |29] have shown that the Tsallis entropy indeed satisfies this criterion 



and hence meets this requirement of concavity. The above features of Tsallis entropy thus 
make it unique among other forms for entropy suggested in the literature. In this context, 
it is very important to understand more deeply the properties of the generalized statistical 
mechanics based on the Tsallis entropy. In particular, a generalization of the approximate 
methods of calculation of its thermodynamical functions is of great value, as for instance, 
the semiclassical approximation |J{J, perturbation and variational methods. The present 
work deals with the last two questions. We develop here the perturbation and variational 
methods in a unified way. This approach provides the generalization of the Feynman proof 
|3T| of the Bogoliubov inequality, which appears to be a natural generalization of this in- 



equality. Indeed, we shall prove that the original form of the inequality is preserved (see 
inequality ([|)). This inequality does not coincide from that proposed in ref. []32]| , except for 



q = 1. This is due to the fact that we use different mathematical inequalities to derive our 
final results. 



The Tsallis entropy and a g-expectation value for an observable A is defined respec- 
tively as S q = kTr p (1 — p q ~ l ) /{q — 1) and (A) q = Trp q A, where p is the density matrix, 
q G R gives the degree of nonextensivity, k is a positive constant. Without loss of generality, 
we employ k = 1 in the following analysis. By using the above definitions with Trp = 1 one 



obtains the canonical distribution [03 



p n = p{E n ) = , (1) 

where {p n } are the probabilities, (3 is the inverse of the temperature, {E n } is the set of 
eigenvalues of the Hamiltonian, and 

Z g = J2[l-(l- q )PE n f^ (2) 

n 

is the generalized partition function. In the expressions (||) and (@) we assumed that 1 — 
(1 — q) f3E n > 0. When this condition is not satisfied we have a cut-off. For instance, when 
a classical partition function is calculated, the integration limits in phase space are given by 
the condition 1 — (1 — q) j3H > 0. ^From Eqs. ([J) and (H) several relations can be obtained, 
for instance, the generalized free energy becomes 

1 z l - q - 1 

F, = U q -TS, = --^—, (3) 

where U q = J2nPnF n is the generalized internal energy, and T is the temperature which 
satisfies the relation 1/T = dS q /dU q . Note that the previous expressions are reduced to the 
usual ones in the limit case q — > 1. 

To develop the perturbation method in this generalized statistical mechanics we assume 
that the Hamiltonian of the system is 

H = H + \Hj. (4) 

In this expression, H Q is the Hamiltonian of a soluble model, XHj is small enough so that 
it can be considered as a perturbation on Ho (Hq and Hi need not necessarily commute), 
and A is the perturbation parameter. Thus, the perturbative expansion of the free energy 
can be written as F q (X) = + XF® + ^-F® + ... . 

To understand how the corrections F^ are calculated it suffices to evaluate the first three 
terms. The first term is the free energy for the case without perturbation, F^ ' = F q (0). 
The second term is obtained from the first derivative of F q (X) at A = 0, 



9F q (0) 
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£[l-(l- g )/3£ r 



,1/(1-9) 



(^)S 0) , 



(5) 



A=0 



<9A dA 

where the superscript (0) indicates that the g-expected value is calculated for A = 0. To 
obtain the last equality it is necessary to exchange the order of derivative with respect to A 
and the sum over n. Furthermore, the Hellmann-Feymann theorem was used , dE n /d\ = 
(n\dH/dX\n), or equivalently the relation d \ n)/d\ = J2 m ^ n ( m \Hi\ n) / (E n — E m ), which 
was obtained from the first order perturbation theory. The calculation of is similar. 
Thus, by using the above considerations we obtain 
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(2) = d 2 F q (0) 
q ' d\ 2 



-fo KT'EM^) (p(E^(n\Hj\n)^ - (Hj) 



S 12 |(n|i?i|m) 



i0) 2 P (E^y-p(E^y 



When Ho and Hj commute (in the classical case, for instance) the expression ([]) is more 
simple, i. e. the second term in the right side is zero. Notice that the Eqs. (|5|) and (|6]) are 
correct for any A, but in this case it is necessary to consider the dependence of | n) with A 
and to substitute E^ by E n . This property will be used later on. 
Finally, the free energy up to the corrections calculated above is 

A 2 



(6) 



Fq W =F q {Q) + \{H I )f-^-(5q{zf)) q - 5>(E<°>) (piE^y-'inmn)^ - (Hj) 



£ \(n\H T \m) 



M 0) 



B?(°) 



+ A 



(7) 



where | n) is evaluated with A = 0. It is important to emphasize that the previous calculation 
was performed supposing that the derivative with respect to A can commute with the sum 
over the states, but generally it is not true because the sum, through its upper limit, can 
depend on A. Indeed, when the factor 1 — (1 — q) ftE n is negative the sum in(Q) must be 
truncated, and as E n depends on A it follows that the upper limit of the sum over the 
states also depends on A. Without losing generality one can consider ft > in the following 
analysis (for ft < the analysis is essentially the same). Therefore, we can suppose that 



E n > for all n. In this case, 1 — (1 — q) (3E n is not negative for q > 1, then the order of 
the derivative in A with the sum in n can be interchanged. For the case q < 1 the analysis 
is more complicated, and we will return to this question later in the classical context. 

Consider now the variational method. In this case, the Hamiltonian of the system is 
H = H + Hi, but in the following analysis it is convenient to employ a Hamiltonian that 
interpolates continuously Ho and H. Thus, we use the Hamiltonian ([|) with A G [0, 1]. Now, 
we will employ the general identity for free energy, 

F q (A) = F q (0) + \F' q (0) + y F" q (A ) , (8) 

where each prime indicates a derivative with respect to A, and Ao G [0,1]. On the other 
hand, the expression (|6|) is not positive for q > 0, and it is not negative for q < for A = A . 
In fact, E„p(#n) [(piE^inlHjln) - (Hj),) 2 ] > 0, and p(E n ) <p{E m ) for E n > E m . 
Therefore, by considering (g) and (|8]) we conclude that F q (A) > F q (0) + \ (Hj) q ^ for q < 0, 
and F q (A) < F q (0) + \(Hj) q ^ for q > 0. But these conclusions, as we discussed at the 
end of the last paragraph, cannot always be true for all q. Thus, when we consider f3 > 0, 
E n > and the fact that the last inequality is correct for all possible values of A, we can 
choose A = 1 and conclude that 

F q < Ff + (Hj)f> (9) 

for q > 1 . As we can see immediately, the inequality (Q) is a natural generalization of the 
usual Bogoliubov inequality, F<F^ + (Hi)(°\ 

To analyse the contribution of the cut-off in the perturbation and variational methods 
it is convenient to consider the generalized classical statistics for (3 > and E n > (the 
analysis for (3 < and E n < can be extended analogously). In this case, there is cut-off 
only for the case q < 1, so we are going to restrict the following analysis to the case q < 1. 
Now, we change J2 n by J dT, where dT = Y\ s (dx s dp s /h) (h is the Planck constant), and 
the integration is over the phase space region defined by the inequality 1 — (1 — q) (3H > 0. 
Thus, to analyse F^ we make use of the following identity 
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dX 



In the present application, / is a function of phase space variables, y u , and dV is the 
hypersurface defined by the equation 1 — (1 — g) (3H = 0. 

When the last term in the right side of (|H]) is zero, the previous analysis is recovered. 
In general, this occurs when / = and dy n /d\ is finite on dV. By using these conditions 
and / = [1 — (1 — q) fiH] 1 ^ 1-9 ^ we conclude immediately that (|5|) remains intact for q < 1. 
For the function / contains a term proportional to [1 — (1 — q) (3H] q ^ 1 ^ q K Thus, the 
condition / = on dV is satisfied only for qj (1 — q) > 0, i. e. g > 0. For F^ the 
function / contains terms proportional to [1 — (1 — q) (]}j'^ nq ~^ n ^ 1 ^^ l ~ q \ Consequently, it is 
necessary q > 1 — 1/n in order that F^ does not contribute to the last term of (|1(]) . These 
observations indicate that arbitrarily high orders in the perturbation expansion, developed 
in preceding discussions, can be used only for q > 1. On the other hand, the variational 



approach can be used, in general , for all q > 0. When the last term of (|10D is not zero, 
there is no guarantee that F^ remains negative, because, in general, the sign of dy u /d\ is 
not known. In this case, the inequality @ cannot be generally employed. 

To illustrate the methods developed above, we consider a one-dimensional harmonic 
oscillator, H = p 2 /2m + mw 2 x 2 /2, which we will approximate by a particle in a square well 
potential, Hq = p 2 /2m + Vq with Vq = for \x\ < L/2 and Vq = oo for \x\ > L/2. In this 
example, we are assuming /3 > 0. The partition function of the unperturbed system and the 
g-expectation value for H\ = H — Hq can be directly calculated, i. e. 
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5 , g< 1 



(11) 



1 >..■'-'•• 2 Mj-fjj q>l 



Thus, the free energy, F 9 (0), can be obtained directly from Eqs. (|TTD and (0). Note that, 
■Z^ ) is convergent only for g < 3. 



Now, we consider the variational method. In this case, the minimization of + (H — 
H )W leads to 

A / q X 1 / 2 

(13) 



3 — q \muj 2 (3 J 

for both q < 1 and q > 1. Substitution of this result in + (H — #o)g leads to the 
optimum approximation for the free energy. Finally, the comparison of this approximation 
for the free energy with the exact one, (valid only for q < 2) 



1 



2tt 1 



- 1 



(14) 



2 - qj 

gives a good approximation, as in the case q — > l(see fig. ([I]) ). Furthermore, fig. ([!]) shows 
that the approximation is improved for larger g. Notice that all previous expressions reduce 
to the usual one in the limit q — > 1. The perturbative contributions can be obtained in a 
similar way. 

Summing up, we have developed here generalized perturbation and variational methods 
for the nonextensive context in a unified way. In this approach, a generalization of the Bo- 
goliubov inequality which is form-invariant for all q is obtained. This property is in variance 



with the generalization presented in ref. [32||. When we consider arbitrary high orders in the 
perturbation expansion, we must have q > 1. On the other hand, the Bogoliubov inequality 
can be generally used for q > 0. We believe that approaches presented here are useful in the 
discussion of the anomalies currently associated with nonextensive systems. 
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FIGURES 

FIG. 1. Free energy approximated and exact vs. temperature for three typical values of q (with 
to = m = h = 1). 



11 



